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The resistive magnetohydrodynamic (MHD) equations as usually defined in the quasineu- 
tral approximation refer to a system of 14 scalar equations in 14 scalar variables, hence are 
determined to be complete and soluble. These equations are a combination of Navier-Stokes 
and a subset of Maxwell's. However, one of the vector equations is actually an identity 
when viewed from the potential formulation of electrodynamics, hence does not determine 
any degrees of freedom. Only by reinstating Gauss's law does the system of equations 
become closed, allowing for the determination of both the current and mass fiow velocity 
from the equations of motion. Results of a typical analysis of the proposed electromagnetic 
hydrodynamic model including the magnetization force are presented. 

Many authors define the low frequency resistive MHD equations 

as the zeroth and first order moments of the Vlasov equation with adiabatic closure in 
conjunction with the two curl equations among Maxwell's. For the neutral fiuid, the sum 
and difference of the ion and electron equations of motion give the net force balance equation 
and the generalized Ohm's law. Using D / Dt = d/dt + Vf ■ V, we write the usual equations: 

^ + V.(p„V,)^0, ^^0, (1) 
Pm^^ = JxB-Vp, r7J = E + V/xB, (2) 

VxE = -^, VxB = /ioJ, (3) 

where r] is the resistivity and 7 is the appropriate index for the case under consideration, 
and the degrees of freedom are pressure p, mass density pm, fiow velocity V/, current J, and 
electromagnetic fields E and B, giving a naive counting of 14 scalar equations for 14 scalar 
variables. However, while for decades [o], [l^ the argument has been made that Gauss's 
law may be neglected with impunity, no one within the plasma physics community has 
denied the applicability of the potential formulation of electrodynamics [l])]. The unnamed 
of Maxwell's equations (often called the "no-monopole" equation, V ■ B = 0) is brought into 



play during t 



equation 
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le determination of plasma equilibrium, via solution of the Grad-Shafranov 



13| in toroidal geometry or otherwise, which by the naive counting of above 
would introduce an additional scalar equation, thus over-determining the system, yet is 
commonly known simply to allow for the expression of the magnetic field in terms of the 
vector potential, B = V x A. The reason doing so is valid is because vector identities 
by mathematical definition do not determine any degrees of freedom; they reduce them. 



Inserting that expression into Faraday's law 



14| . we recover V x (E + dA/dt) = 0, whence 



E = — V$ — dA/dt, which clearly displays the division of the electric field into static and 
dynamic components and reduces three of our naive degrees of freedom down to one for which 
we have no equation. Unless one wishes to invent new physics, the resolution is clear — the 
reinstatement of Gauss's law, V ■ E = — — (9(V ■ A)/dt = Pe/^o which vanishes for a 
neutral fiuid, is required to close the system of equations, bringing the number of scalar 
equations and degrees of freedom into agreement with the number 14 — 3 + 1 = 14 — 2 = 12. 
We remark that Faraday's law is no less an identity than the no-monopole equation as both 
are given by the general theory of vector fields. Gauge invariance plays a special role in the 
local conservation of charge, best expressed in manifestly Lorentz covariant notation. 

From a particle physicist's field-theoretic point of view 15, 16, 17, [is], [3], the Maxwell 
field tensor F'^'^ = d^A'^ — in media is known to have only 3 physical degrees of free- 
dom embodied by the four-potential A^ = ($/c. A) subject to the gauge condition, not 
3 for each of the electric and magnetic fields, which couple to sources given by the con- 
served four-current = {cpe, J) through the inhomogeneous Maxwell equations d^F^'^ = 
{dnd^)A'^ — d'^ld^A'^) = fioJ"^ , which are explicitly Lorentz covariant and also gauge invari- 
ant, and the homogeneous Maxwell equations, given by the divergence of the dual tensor 
pfj.u ^ e^""^ Fai3/2, where e^'^"^ is the permutation tensor, as d^F'^'^ = 0, are satisfied iden- 
tically when written in terms of the electromagnetic potential, hence do not determine any 
degrees of freedom. Antisymmetry in F'^'^ immediately implies conservation of the current, 
dydfj^F^^ = poduJ'^ = 0, thus it carries only 3 degrees of freedom also. One may recast the 
Maxwell equations into a component-free form through the use of differential geometry 18| , 
where "the existence of integrals implies a duality between forms and chains" which may be 
exploited. In natural units /ig = cq = c = 1 and using the exterior derivative d, the Hodge 
dual *, the connection 1-form A = A^dx^, the curvature 2-form F = {—F^j_^/2)dx^ A dx'^ , 
and the current 3-form J = {Jxdy A dz + Jydz /\ dx + J^dx A dy) A dt — Pedx A dy A dz 
which satisfies the continuity equation d J = 0, one writes the field equation as d *F = J 
and the Bianchi identity, which is a statement on the structure of the manifold, as d F = 0, 
whence F = dA, and we remark that gauge invariance, through Noether's theorem, implies 
conservation of the covariant current. What all this shows is that the natural, physical divi- 
sion of the Maxwell equations is not into the divergence and curl equations but rather into 
the homogeneous and inhomogeneous equations, whereby the Bianchi identity carries the 
structure for the potential formulation and the field equation carries the dynamics obtained 



3 



from the action. 

The imphcation for plasma physics is clear: the quasineutral approximation does worse 
than just neglect an effect, as it introduces inconsistency into the equations when the com- 



ponents of the electrostatic field are treated in isolation [20|]. Arguing that Maxwell's diver- 
gence equations are initial conditions for the curl equations is incorrect in media, for while 
in vacuum such statement leads to the propagation of electromagnetic radiation with two 
physical states of polarization, the source terms spoil such interpretation, and the divergence 
of the Maxwell-Ampere equation only recovers the equation for local charge conservation, 
which must be respected independently of the conservation of mass addresed by the zeroth 
moment of the Vlasov equation, when Gauss's law retains its intended form. Note that 
authors not including the no-monopole equation explicitly within the system do not make 
the argument of having 14 equations and degrees of freedom, as that equation represents an 
additional member. Claiming that in general the sources may be uniquely determined from 
expressions for the fields is inappropriate, for while suitable boundary conditions must be 
supplied, the differential operators hence the boundary conditions are applied to the fields, 
not the sources. The reason for the expression "Maxwell-Lorentz electrodynamics" is be- 
cause Maxwell's theory tells one how the fields react to the sources, and the Lorentz force 
through the equations of motion tells the sources how to react to the fields; trying to go the 
other way around the loop is not well defined, as the physics is contained within the action 
from which both the field and source equations of motion may be obtained. 

Let us examine in detail where difficulties are encountered by the neoclassical approach, 
a term we use to encompass all non-classical approaches to the fluid description of ion- 
ized particles regardless of geometry — such discussion [2l| invariably engenders a hostile 



response 
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231] from its adherents yet is necessary if one is to consider the application of 



electro dynamic field theory in tensor notation to the many-body system commonly called 
a plasma. The scalar degrees of freedom pm and p may be associated with the scalar equa- 
tions for mass and energy conservation. Equations ([T]), as no other quantities appear in 
those equations for the case of vanishing flow velocity; the presence of a flow velocity Vj 
couples those equations to the rest of the system to be solved simultaneously. Note that the 
previous argument of the second paragraph tacitly assumed that the equations of motion 
in the form of the generalized Ohm's law and the convective force balance. Equations ([2]), 
were associated with the degrees of freedom {V/, J}; whereas here, without Gauss's law, one 
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must determine the electric field from a n eq uation of motion, usually the generalized Ohm's 
law (however the ion 2J] and electron 25|] equations of motion are also used), giving the 
solution E„eo = rjj — Vf xB. Faraday's law in conjunction with the no-monopole equation 
then relates the electric field to the potentials — E„eo = dA/dt + V$, where without Pois- 
son's equation or its gauge invariant generalization the relation between the potentials and 
the space charge density pe remains unspecified (in essence, Faraday's law here determines 
a potential $ which is not an independent degree of freedom), and its divergence gives in 
various gauges 

d_ 
di 



V ■ (V/ X B - r/J) 



-V- A + 



Coloumb (V • A = 0) 

d 

Lorenz (V ■ A = -pQtQ—^) 
Weyl ($ = 0) 



□2$ 
d 



dt 



V- A 



(4) 
(5) 
(6) 

(7) 



where the LHS is explicitly gauge invariant whereas the form and interpretation of the RHS 
is dependent upon one's choice of gauge. The issue of gauge invariance is a red herring in 
the discussion, for while true physics must be equally described in any and all gauges, the 
crucial error in the neoclassical approach is its use of an equation of motion to determine 
the electric field, which does not respect Lorentz covariance. (Note that modern power 
generators and electric motors certainly are not moving materially at relativistic speeds yet 
make full and practical use of the covariant transformation properties of the field tensor 
through Faraday's law of induction.) Returning to the expression for E„eo, let us now 
examine its transformation properties under a change of reference frame. Let 5* be the 
frame of the neoclassical observer, and let S' be the frame moving with velocity Vj with 
respect to S. Without loss of generality, the flow velocity in S is taken along the x-axis, thus 
V/ = {Vf, 0, 0) =^ gives E„eo = {il^x.^Jy + ^fSz^i^Jz — VfBy), using Einstein's velocity 



addition rule 26| gives = 0, and for 7 = 1/ yl — V'j jc? the transformation for proper 
velocity applies to J, the spatial part of the four-current J^^^ = (0, J), giving 









VJx 


Keo = VJ' = 






IVJy 











E' 



(8) 



where E' = [Ex,'j{Ey — VfBz),'j{Ez + VjBy)] is the transformation law for the physical 



5 



electric field. Equality could hold only if r^' = 77/7 = 777 implying 7 = 1, which holds only 
when Vf = 0, thus only in the neoclassical frame of reference but also implying a vanishing 
flow velocity, contradicting the initial assumption V/ 7^ 0. The expression for E„eo has 
inherited the nature of a velocity vector from its neoclassical determination hence cannot 
possibl y re present a true electric field, which does not transform as the spatial part of a four- 



vector [11, 14, 26|. Furthermore, as ultimately B(J) may be determined from Ampere's law 



or the equivalent Biot-Savart law (for steady currents only as all the terms with E need be 
present for Maxwell's theory to respect local charge conservation), the neoclassical electric 
field depends explicitly on the two vectorial quantities of current and mass flow, E„eo(J, V/). 
In order to completely determine the system, both of those quantities must find solution; 
however, having already used one of our equations of motion in the guise of Ohm's law, we 
have left only one vector equation for the net conservation of momentum, pmDV f / Dt+Vp = 
J X B, which leaves one vector's worth of degrees of freedom without solution, leading to the 
use of a stationary equilibrium equation Vp = J x B in the analysis of non-stationary plasma 
experiments 24, l25|. We note that the predictions of the neoclassical (NCLASS) model for 
the poloidal velocity found in a tokamak presented in Reference I2J] explicitly fail to agree 
with the experimental measurements. By reinstating the determination of the electrostatic 
field via Gauss's law, what returns is the generalized Ohm's law, an equation of motion 
which one may solve for the motion appearing in that equation, which in conjunction with 
the convective force balance equation fully determines the system. Ultimately, the various 
arguments presented in support of the neglect of Gauss's law are superseded by the rigorous 
formalism of differential geometry, whereby casting the Maxwell equations into intrinsic, 
geometric form, d*dA = J, comprises very deep and powerful statements concerning what 
is known about our Universe. 

Consequently, we here advocate the use of the electromagnetic hydrodynamic (EMHD) 
equations specified by the fluid equations for all species s, 

^""^""^ + V ■ (n.m.V.) = , f I- + V. ■ = , (9) 



dt y ^ ^ -J ' \dt ^ J {nsmsp 

risTTis + ■ V j V, + Vps = n,e, (E + V, x B) + ^ F.^ , (10) 

^ ' k 

associated with the degrees of freedom {ns,Ts, V^} for ps = nJFs = nsksTs, where k ^ s 
for the friction term F^fc = — F^^ representing interspecies collisions and neglecting viscos- 
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ity, polarization, and magnetization, in conjunction with the field equation d^F^'^ = j^oJ" 
associated with the degrees of freedom A^. Extension to incorporate the magnetic decom- 
position B//io = H + M into free and bound currents is straightforward, and a gyrotropic 
pressure tensor is forthcoming, noting that the fiuid mechanics of Equations and flTU]) 
remains to be cast into the Lorentz covariant form manifest for the electromagnetic sector. 
As a preliminary, we solve the Equations (ITOl) for everyone's favorite hydrogenic plasma of 
fusion interest, fully ionized deuterium, in the neutral fiuid limit pe — > such that V-E ^ 0, 
for the case of an infinite column at equilibrium d/dt ^ with prescribed coaxial applied 
electric and magnetic fields and an assumed pressure profile. With the definitions for species 
s G {e, i} of particle density n = rie+rii = 2no, mass density pm = Xls ''^s^s, free momentum 
density PmYf = J^s^s^sYs, and free current density J/ = Xls'^'s^sV^, one finds 



Ve 



epn 



-rrii 
rrie 



(11) 



where Vg^j are the species' fiuid or guiding center velocities. Gyromotion is so far incor- 
porated through the inclusion of the drift momentum pmVd = V x as the curl of the 
angular momentum density carried by the cyclotron motion and the diamagnetic current 
= V X M as the curl of the magnetic moment per unit volume, giving total momentum 
density pmV = pm(V/ + V^) and total current density J = J/ + J^, and through the mag- 
netization force Fa/ = poVM • H = -p^oVMH as M = -Mh = -(M/i7)H, an important 



effect neglected in the usual analysis o: 
tions in fusion 
science 
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371, 
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magnetic fiuids 
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ary equilibrium elsewhere ■ 
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alasma equilibrium despite its exp erimental applica- 
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471] and addressed for the case of a station- 
. For the gyro-momentum, we take Lg = Lgg = J2s 



for Is = msv'j_gUJs/ujg = 2Tsujs/ojl where Us = — e^Bs/m^ and v±s = \/'2W^ jnis^ and for the 
species dipole density = n^ps = —nsiyV^jBi^'h, where the total field felt by a particle of 
species s is B^/po = H + M — p^ = H + Mfc + asM^ ^ {H — M)h for = (n^ — l)/?2s ^ 1 at 
sufficient plasma density and points along bg = Bs/i?s = h, giving a total magnetic dipole 
moment of M = J2s fs^s = 2no(/ePe + /i/ii)/2 = —njlh. The coUisionality factor /, defined 
using 

/. = e^(l-e-/-)|"' for..,».., 

I -> uJs/ Us for < z/^ , 



(12) 
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where Us = Vss + ^sfc, represents collisional disruption of the gyromotion and is normahzed 
so that fp = J2s fsPs- Using p = p/ fiQ, the general magnetization model 

fePe , fiPi fp 



M 



H - Mi- aoMe ^ H - Me- a^Mi H - aM 



(13) 



has solution M/H = (1 - y^l - iafp/ H^)/2a for < a < 1 and /5-limit P = 2p/H^ < 
l/2af on the ratio of the kinetic to the magnetic pressure, which goes to 1/2 for a, / — 1 
and to 1// for a —>■ 1/2. As uJs is in terms of Bg, an iterative approach to the collisionality 
factor may be defined by Mn{H — aMn) = fnP, starting at /o = 1 with /„ = /„(M„_i), and 
for a dense, magnetically confined plasma we find / remains very close to unity. 
The sum of Equations (fTO!) gives the equilibrium net force balance equation 



C+ + Vp+ = J X B + /ioV(M+ ■ H) 



now including Fm, and their difference the generalized Ohm's law equation 



C_ + Vp_ = rioe (Vi + V^) x B + /xoV(M_ ■ H) + 2 [riocE - 



(14) 



(15) 



where noe(Vj + Vg) = 2noeV — no(mj — me)J/pm and {Vf,3f} are replaced by {V, J} 
everywhere except in the friction term, p± = pi ± pe and M-t = /^Mj ± /eMe, and the 
convective terms C+,_ = tiq [nii (Vj ■ V) ± rrie (Vg ■ V) Vg] are given by 



C_ 
C. 



no {rrii + rrie) 

no {iTLi - rrie) 
2nomemi 



V ■ V V + — ■ V 



(V ■ V) V 



mprrij 



■ V 



Pr. 



Pm, 

3_ 

Pm, 



[V ■ V) — + f — ■ V ) V 

Pm V Pm 



(16) 
(17) 
(18) 



Specializing consideration to an infinite, axisymmetric plasma column in (r, 6, z) coordinates 
with V — (9/9r, 0,0) and applied fields Eq = Eqz and Bq = Bqz, we define profiles 
for the density and temperature with profile parameter ax as no(r) = no(0) + [no(ri) — 
no(0)](r/ri)'^"o for ri at the boundary of consideration and similarly for T^a- Equilibrium 
requires divergence- free momentum and current, V ■ (pm^) = V ■ J = 0, and the convective 
terms reduce to (A • V)B — > — AeSgr/r^. The 9 and z components of Equation fllSI) give 
the free current in terms of the applied electric field J/ = {Q^Q^nQe^EQ/rrieVei)-, from which 
the enclosed current Iz{,f) = dr'2TTr'Jfz{r') determines the free field Hf{r) = Iz{r)/2TTr 
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such that H = Hq + HfO. Our logic follows Eq — J/ — H/ + Hq — M — > J(i, where the 
diamagnetic current is given by 

Jrf = V X M = -V X (^^H^ = H X - X H , (19) 

and similarly for 'Vd- Scattering times Tgk = ^/vsk are calculated using common formulas 49|, 



50( 1 ■ The mass flow velocity V is then found from the radial components of Equations ([T 



and (fTSll . Results for a typical analysis are found in Figures (1) and (2) with the numerical 
parameters given in the first figure caption. Note that the theory above accounts naturally 
for the "spontaneous" generation of a flow velocity as a consequence of the equilibrium force 
balance equations and offers an alternative to the neoclassical approach for the analysis of 



plasma experiments 



24 



25|. 



In conclusion, we have found that Gauss's law may not be extracted from the covariant 
field tensor equation without introducing serious mathematical inconsistencies. The expres- 
sion of Maxwell's theory in intrinsic, geometric form holds at the classical level, and suitably 
quantized its predictions have been verified to an accuracy of parts per billion. To conjecture 
a range of energy density at which it fails to hold is inconsistent with experience of both the 
everyday world and precision experimental tests. Reinstating Gauss's law lets one predict 
both the momentum and current within a plasma from the principles of electromagnetic hy- 
drodynamic theory, allowing for the experimental testing of the hypothesis that Maxwell's 
equations really do apply to the physics of many-body systems. 
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FIG. 1: (Color online.) Results from a typical analysis of magnetically confined, fully ionized 
deuterium with ri = Im, applied electric field £"0 = 20mV/m, and applied magnetic field Bq = 2T. 
(a) Input pressure profiles with central density no(0) = 9 x 10^^ /m^, central temperature T(,^i{0) = 
3keV, profile parameter ax = 2, and pedestal parameter X{0)/X{1) = 10. (b) Calculated /3-limit 
and (3 with a central value of 5.4%. (c) The magnetization force equals or exceeds the pressure 
gradient force in the equilibrium equation, (d) The net magnetic field is reduced from the free 
current value by the magnetization. 
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FIG. 2: (Color online.) Results from a typical analysis of magnetically confined, fully ionized 
deuterium with ri = Im and total enclosed free current = 2.8MA. (a) The diamagnetic current 
has both an azimuthal and axial component, (b) The net axial current is reduced from the free 
current value by the axial diamagnetic current, (c) & (d) The effect of the gyromotion on the plasma 
flow velocity is small but non- negligible. Here, Vj is the flow with neglect of the gyromotion and 
V includes the drift velocity V^, found to be on the order of mm/s. 
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